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Abstract 

Two nozzle performance prediction procedures 
whien are based on the standardized JANNAF meth- 
odology are presented and compared for four 
rocket engine nozzles. The first procedure 
required operator intercedance to transfer data 
between the individual performance programs. The 
second procedure is more automated in that all 
necessary programs are collected into a single 
computer code, thereby eliminating the need for 
data reformatting. Results from both procedures 
show similar trends but quantitative differences. 
Agreement was best in the predictions of specific 
impulse and local skin friction coefficient. 

Other compared quantities include characteristic 
velocity, thrust coefficient, thrust decrement, 
boundary layer displacement thickness, momentum 
thickness, and heat loss rate to the wall. 

Effects of wall temperature profile used as an 
input to the programs was investigated by running 
three wall temperature profiles. It was found 
that this change greatly affected the boundary 
layer displacement thickness and heat loss to the 
wall. The other quantities, however, were not 
drastically affected by the wall temperature pro- 
file change. 

Introduction 

Accurate predictions of boundary layer 
growth are necessary for accurate predictions in 
rocket engine performance. Equilibrium and chem- 
ical kinetics programs already exist which pre- 
dict performance losses with sufficient accuracy. 
However, the computer modeling of boundary layer 
growth and its effect on rocket engine perform- 
ance is less standardized. This is due to dif- 
ferences in the physical modeling and numerical 
procedures used. Boundary layer losses start to 
detract greatly from the overall performance in 
rocket engine nozzle designs with large exit-area 
to throat-area ratios. These high area ratio 
nozzles (>400) are of interest for orbital trans- 
fer vehicles to obtain maximum performance as 
indicated by high specific impulse. 

The purpose of this paper is to compare the 
results from two procedures for evaluating the 
performance of nozzle contours. Both procedures 
are based on the standardized JANNAF methodology* 
for predicting rocket engine nozzle performance. 
The first, the disjoint procedure, requires 
transmission of data between three separate 
performance prediction programs. The boundary 
layer prediction program in this procedure is 
BLIMP. 2 BLIMP is an independent, stand-alone 
program and can only be utilized after prior runs 
of other performance programs. In this case, two 
programs were run as an input to BLIMP. The 
second, the integrated procedure, utilizes one 
performance program which has many subroutines to 
perform the prediction. BLM 3 was the boundary 


layer prediction program used in this procedure. 
Unlike BLIMP, BLM is accompanied by a set of sub- 
programs which will calculate all the necessary 
preliminary information required for the boundary 
layer analysis. This total set of subprograms 
(including BLM) is referred to as the 1985 ver- 
sion of the Two-Dimensional Kinetics program 
(TDK. 85). 4 

Performance comparisons were made for noz- 
zles having area ratios 60:1, 200:1, 400:1, and 
1000:1. Each nozzle had a throat diameter of 
1 in. and was specified to run at a chamber pres- 
sure of 1000 psia using hydrogen and oxygen as 
the propellants. Variations of the wall temper- 
ature profile were made to investigate its effect 
on the engine performance. 

Description of Procedures 

The first procedure for evaluating the 
rocket engine nozzle performance is referred to 
as the disjoint procedure. The succession of 
programs used in this procedure is the One- 
Dimensional Equilibrium (ODE) program, 5 the 
One-Dimensional Kinetics (ODK) program, 5 and the 
Boundary Layer Integral Matrix Procedure program 
(BLIMP). 2 ODE is a thermochemical oroperties 
program which provided throat chemical composi- 
tion input data for the ODK program. ODK, a 
kinetics program which computes finite rate reac- 
tion performance, was run from the nozzle throat 
to the exit. The results provide BLIMP with pres- 
sure profile and chemistry information as input. 
The exclusion of a two-dimensional kinetics anal- 
ysis results from the fact that the 1973 version 
of the Two-Dimensional Kinetics program (TDK. 73) 
was to be used, but it gave unstable results. 
Therefore, an average value for the equivalent 
two-dimensional divergence efficiency of 0.996 
was assumed in lieu of running TDK. 73. This 
divergence efficiency is equivalent to a 7.25° 
divergence half-angle.' The flow diagram for 
this procedure is shown in Fig. 1. Table 1 shows 
the required inputs for each program used in the 
disjoint procedure and from which program (if 
any) they are generated. 

The second procedure used is called the inte- 
grated procedure. In this study, the integrated 
procedure refers to a single program, the 1985 
version of the TDK program (TDK. 85). TDK. 85 com- 
bines TDK. 73 as a subprogram along with the 
boundary layer module (BLM), ODE, ODK, and many 
other subprograms into one large program. The 
inputs required for the integrated procedure are 
shown in Table 2. 

Description of Programs 

The individual programs described below are 
performance prediction codes. 
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One-Dimensional Equilibrium Program (ODE) 5 

This program calculates the equilibrium con- 
ditions for a reactive mixture in an idealized 
closed combustion system. Reaction rates are 
considered to be infinitely fast so that all reac- 
tions reach local equilibrium conditions. A 
rocket performance option is run which predicts 
the theoretical maximum performance of isentropic 
convergent-subsonic and di vergent-supersonic flow 
in a nozzle with shifting equilibrium chemical 
composition during expansion. The results from 
ODE were used as throat condition inputs to the 
ODK program in the disjoint procedure. 

One-Dimensional Kinetics Program (ODK) 6 

This program calculates the deviation from 
equilibrium performance in the nozzle due to 
finite reaction rates not allowing chemical equi- 
librium to be reached. Unlike ODE, ODK bases the 
fraction of mixture which will react on the reac- 
tion rates of the mixture. ODK considers complex, 
homogeneous, ideal gas reactions in combustion 
and nozzle expansion situations. Eighteen chemi- 
cal reactions and four three-body recombination- 
dissociation reactions were used by ODK for the 
fuel /oxygen combination of hydrogen and oxygen 
(Table 3). The output from this program provides 
an axial pressure profile along a streamline 
boundary (i.e. the nozzle wall) as an input to 
the BLIMP program in the disjoint procedure. ODK 
also outputs performance results including speci- 
fic impulse, thrust coefficient (Cp), and char- 
acteristic velocity (c*). 

Boundary Layer Integral Matrix Procedure Program 

This program calculates the boundary layer 
growth in a rocket nozzle. For axi symmetric 
flow, the program can calculate either shifting 
equilibrium or frozen compositions for any fuel/ 
oxidizer combination. The program uses an inte- 
gral matrix solution algorithm which is equivalent 
to a higher order finite difference approach. 

This program will output various boundary layer 
properties as well as the coordinates of the 
potential flowfield if the input described the 
real wall, or the coordinates of the real wall if 
the potential flow contour was input. 

Two-Dimensional Kinetics Program 1983 Version 

rrDK'.SS ) 3 

This program is really a combination of pro- 
grams. It incorporates the ODE and ODK programs 
along with a transsonic module (TRANS), a method 
of characteristics routine (MOC), a two- 
dimensional kinetics effects module, and a bound- 
ary layer predictor (BLM). TDK. 83 integrates all 
the pieces which were used in the disjoint proce- 
dure and passes the needed information between 
the programs automatically. The results of each 
run are printed out in a summary table at the end 
of a run. The program provides the user with the 
capability of using specific combinations of sub- 
routines for each section of the engine. All the 
computations involved with analyzing a rocket 
engine nozzle are achieved with one computer run. 
This includes evaluation of the kinetic devia- 
tions from equilibrium performance, calculation 
of the sonic line in the throat region, setting 
up the mesh of characteristic lines, calculating 


the kinetics variances in the nozzle downstream 
of the throat, and predicting the boundary layer 
growth and losses. Also, the program can dis- 
place a potential wall contour by the displace- 
ment thickness calculated by BLM and recompute 
the kinetic variances in the nozzle to attain a 
better estimate of the boundary layer thickness. 

Description of Input Cases 

Four nozzles, with area ratios of 60:1, 

200:1, 400:1, and 1000:1, were chosen for this 
evaluation study. All four had 0.5 in. throat 
radii and lengths corresponding to a 100 percent 
15° cone with the same throat radius and area 
ratio. ' They all operated at 1000 psia chamber 
pressure and used normal boiling point (NBP) liq- 
uid hydrogen and NBP liquid oxygen as their fuel 
and oxidizer. The propellants were injected with 
a mixture ratio (0/F) of six. 

For the given area ratios, the nozzle exit 
radii were determined. These exit plane lengths 
and radii values are tabulated in Table 4. Since 
the purpose of this study was comparative in 
nature, the actual contour of the nozzles was not 
important so long as they were the same for both 
procedures. The contour for each nozzle between 
the throat and exit plane was calculated using 
the RAO method, a method of characteristics rou- 
tine coupled with an optimization scheme based on 
the calculus of variations. ^ These contours 
(Appendix A) defined the potential flow bound- 
aries for each nozzle. 

The remaining operating conditions of the 
engine were used as inputs for ODK, BLIMP, and 
TDK. The energy level associated with the oxygen 
is -2948 cal/mole and -1977 cal/mole for the 
hydrogen. The geometry variables of the combus- 
tion chamber, throat, and nozzle are shown in 
Fig. 2 and the corresponding values are specified 
in Table 5. 

For both evaluation procedures, the poten- 
tial wall contour was defined by specifying axial 
location values and radii at equal intervals 
along the nozzle centerline. In the integrated 
procedure, the number of contour points was 
greatly reduced at the higher area ratios (>200:1) 
to insure that the contour interpolating routine 
would not compute an erratic wall slope between 
input wall points. The number of points given 
for each area ratio for both procedures, the num- 
ber of points on the sonic line at the throat 
(required to set-up the method of characteristics 
mesh), and the number of segments used in the 
boundary layer prediction (for the integrated 
procedure) is shown in Table 6. 

Each nozzle had the same baseline wall tem- 
perature profile. The wall temperature profile 
of the 1000:1 nozzle is shown in Fig. 3. The 
smaller area ratio nozzles used this same profile 
truncated at the axial position of the nozzle 
exit. The profile for all four nozzles assumed a 
cooled throat section, and a radiation cooled 
nozzle. Near the throat, the profile used is 
very similar to that calculated for the Space 
Shuttle Main Engine (SSME). 2 It is much cooler 
in the throat region than the calculated tempera- 
ture for the RLIO^ due to the difference in 
materials of construction. Figure 4 shows the 
wall temperature profiles of the 200:1 nozzle, 
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the RL10, and the SSME superimposed on the same 
plot for comparison purposes. 

Discussion of Results 

The results of the computer runs showed sim- 
ilar trends for all the area ratios considered. 

The general discussion of results will cover all 
the nozzles investigated, however, only figures 
for the 1000:1 case will be shown since they are 
characteristic of all the results. The 1000:1 
area ratio is also of particular interest for an 
orbital transfer vehicle. The quantities which 
are compared include the boundary layer displace- 
ment thickness, momentum thickness, local skin 
friction coefficient, thrust loss, heat loss to 
the wall, the characteri stic velocity (c*), thrust 
coefficient (Cp), and specific impulse ( I S p) - 

Boundary Layer Displacement and Momentum Thickness 

Figure 5 is a comparison between the bound- 
ary layer displacement thickness as calculated by 
the disjoint procedure and the integrated proce- 
dure for the 1000:1 nozzle. The boundary layer 
displacement thickness predictions for the two 
procedures agree near the throat but then diverge 
monotonical ly toward the exit plane. Agreement 
between the two code's exit plane predictions is 
very good for the 60:1 area ratio case (only a 
1.4 percent difference). However, the 200:1, 
400:1, and 1000:1 cases show a 10.5, 39.6, and 
35.8 percent difference between the exit plane 
values for the two procedures, respectively. The 
disjoint procedure predicts a greater displace- 
ment thickness value and rate of growth than does 
the integrated procedure. These differences 
between the two procedures in boundary layer dis- 
placement thickness are good examples of the rea- 
son for this study. 

Figure 6 shows a comparison of the momentum 
thickness profiles in the axial direction as pre- 
dicted by the disjoint procedure and the integra- 
ted procedure for the 1000:1 area ratio nozzle. 
The integrated procedure predicts a nearly linear 
increase in momentum thickness from the throat to 
the exit, while the disjoint procedure shows a 
gently concave profile which is always greater 
than the integrated procedure profile. This 
results in a slightly higher prediction by the 
disjoint procedure at the exit plane for every 
case. Differences between the two procedures' 
momentum thickness predictions are 11.6, 11.3, 
18.4, and 16.2 percent for the 60:1, 200:1, 

400:1, and 1000:1 area ratio cases respectively. 

A possible reason for the discrepancies 
between the two procedures is that the wall 
pressure profile resulting from the disjoint pro- 
cedure was different than that used in the inte- 
grated procedure. In the disjoint procedure, the 
wall pressure profile was output from the ODK 
program and used as an input to the BLIMP program. 
In the integrated procedure, the wall pressures 
were calculated internally by the TDK subprogram. 
TDK will generally predict larger values of wall 
pressure than will ODK. The differences between 
these values for similar axial stations were as 
large as a factor of two. This large discrepancy 
can change the fluid properties at the edge of 
the boundary layer drastically. This, in turn, 
affects the boundary layer profiles as seen in 
the displacement or momentum thickness plots. 


More subtle differences are found in the 
eddy-viscosity models used by the two boundary 
layer programs. These also directly affect the 
prediction of the boundary layer profiles. BLIMP 
uses the original Cebeci-Smith eddy-viscosity 
model whereas BLM uses a revised version of the 
same model. For example, BLM uses kinematic vis- 
cosity, intermittancy, and transitional terms in 
the near-wall region expression of the dimension- 
less eddy-viscosity whereas BLIMP does not. 

Local Skin Friction Coefficient 

Figure 7 shows the local skin friction coef- 
ficient for the 1000:1 nozzle. The local skin 
friction coefficient is defined as the ratio of 
shear stress at the wall to the local dynamic 
pressure. The closest agreement at the exit plane 
(a difference of 1.5 percent) occurs in the 1000:1 
area ratio case and the largest difference 
(10.2 percent) occurs at the 400:1 area ratio. 

Thrust Decrement Due to Boundary Layer Growth 

The thrust decrement is the loss of thrust 
due to skin friction effects on the boundary layer 
displacement thickness and the wall pressure pro- 
file in the nozzle. Figure 8 shows the predic- 
tion comparisons of thrust decrement for the 
1000:1 case. In all four area ratio cases the 
prediction at the throat by the disjoint proce- 
dure is greater than that of the integrated pro- 
cedure. The integrated procedure diverges from 
the disjoint procedure solution downstream of the 
throat and predicts a greater value at the exit 
plane, as is shown in the figure. The thrust 
decrement predictions vary by 12.5, 20.3, 22.3, 
and 32.2 percent for the 60:1, 200:1, 400:1 and 
1000:1 nozzles, respectively. 

Since the thrust decrement can be directly 
related to the size of the boundary layer, the 
reasons for the discrepancies between the two 
procedures are similar to those described in the 
displacement thickness section. 

Heat Loss to the Wall 

The heat loss to the wall profiles are 
derived from the integration of the amount of 
heat being directed at the nozzle wall per unit 
area per unit time. By integrating this quantity 
locally over the entire nozzle length, the result 
is a cumulative heating rate, with units of 
8TU/sec. The value of heat loss to the wall is 
primarily dependent on the wall temperature pro- 
file and the heat transfer coefficient. Integra- 
tion is performed by the computer programs. The 
comparison of the profiles of the total heat loss 
rate to the wall are characterized by the 1000:1 
area ratio case, as shown in Fig. 9. The values 
predicted by the disjoint procedure are plotted 
starting from the throat whereas the integrated 
procedure's values begin approximately 20 percent 
of the nozzle length downstream of the throat. 

This is because the integrated procedure divides 
the nozzle into segments and gives an integrated 
heat flux value at the end of each segment. The 
values predicted by both procedures at this point 
show good agreement in the throat region. How- 
ever, as the analysis travels down the nozzle 
toward the exit plane, the predicted values 
diverge rapidly. The difference between predic- 
ted values at the exit plane are 14.9, 20.0, 22.7, 
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and 26.8 percent for the 60:1, 200:1, 400:1, and 
1000:1 nozzles respectively. In all cases, the 
integrated procedure predicts a greater value of 
heat loss to the wall for the same wall tempera- 
ture profile. The difference is probably due to 
the difference in wall pressure profiles used in 
the two procedures as well as a difference in the 
surface coefficients of heat transfer. 

Characteristic Velocity, c* 

Analytically, the characteristic velocity, 
c*, is a quantity which reflects the effective 
energy level of the propellants and subsonic 
expansion. Values for the characteristic velocity 
as predicted by both procedures are tabulated in 
Table 8. The characteri stic velocity values pre- 
dicted by the ODE portion of both procedures are 
constant for all area ratios. There is an almost 
negligible (0.2 percent) difference between the 
values calculated by the two procedures, with the 
disjoint procedure predicting a lower value than 
the integrated procedure. Differences between 
the thermodynamic databases used in the disjoint 
and the integrated procedures ' ODE versions can 
account for the difference in these predicted 
values. 

At the boundary layer level of computations 
(accounting for the area change at the throat), 
the difference between the computed c* values are 
within 0.6 percent of each other. Again, the 
integrated procedure predicted a greater value 
than did the disjoint procedure. 

The thermodynamic databases used in the inte- 
grated and disjoint procedures were different. 

This is most likely the cause of the differences 
between the two procedures. 

Thrust Coefficient, Cp 

Tabulated values of the thrust coefficient 
can be found in Table 8. The thrust coefficient 
as calculated by ODE in both procedures show very 
good agreement. The variation in equilibrium 
values between procedures is less than 0.1 percent 
for all area ratios. This variation is caused by 
the thermodynamic database difference mentioned 
in the previous section. 

In all cases, the integrated procedure con- 
sistently predicts a slightly higher Cp value. 

The thrust coefficient for the two different pro- 
cedures differs by 4.7, 3.8, 3.2, and 2.4 percent 
for the 60:1, 200:1, 400:1, and 1000:1 area ratio 
cases, respectively. 

This result taken with respect to the heat 
loss rate results point to an overall consistency 
between the two procedures. The disjoint proce- 
dure predicts a lower heat loss rate as well as a 
lower thrust coefficient. 

Vacuum Specific Impulse, I Q p 

The vacuum specific impulse is a measure of 
the efficiency of a rocket engine. The thrust 
decrement is the loss of thrust due to boundary 
layer growth in the rocket nozzle. These two 
quantities are discussed together because the 
thrust decrement is converted to an impulse decre- 
ment and subtracted from the ideal value of I sp 
predicted prior to the boundary layer analysis. 


The absolute specific impulse values calcu- 
lated by both procedures are all very close, as 
shown in Fig. 10. The predictions differ by 0.10, 
0.51, 0.64, and 1.30 percent for the 60:1, 200:1, 
400:1, and the 1000:1 case, respectively. In all 
cases, the disjoint procedure predicts an I sp 
value higher than those of the integrated proce- 
dure. The specific impulse calculated by the ODE 
portion of both procedures is within 0.25 percent 
due to thermodynamic database differences men- 
tioned previously. The integrated procedure con- 
sistently predicts the larger value. 

The discrepancies in the values of thrust 
decrement predictions translate directly to dis- 
crepancies in specific impulse loss discrepancies. 
The disjoint procedure predicts an I sp loss of 
11.0, 12.4, 13.0, and 13.6 sec from the 60:1 noz- 
zle to the 1000:1 nozzle. This contrasts sharply 
with the 12.7, 15.6, 16.8, and 20.2 sec losses 
which the integrated procedure predicts for the 
60:1 to 1000:1 nozzles. 

The disagreement in specific impulse at the 
1000:1 area ratio result (Table 8) is surprisingly 
low, especially in light of the large thrust 
decrement discrepancy between the two procedures. 
The two orocedures were in good agreement on 
specific impulse predictions. 

Notes on Effects of Wall Temperature Profile on 
Res ults 

In an attempt to investigate the effect of 
wall temperature profile on the predictions, the 
200:1 area ratio case was rerun using the 
integrated procedure with two different wall tem- 
perature profiles. The two new wall temperature 
profiles were the baseline wall temperature pro- 
file displaced by ±500 degrees Rankine (referred 
to as +500 °R and -500 °R wall temperature pro- 
files, respectively). All three profiles are 
plotted in Fig. 11. 

The effect of wall temperature profile on 
the momentum thickness, thrust decrement, and 
local skin friction coefficient is small (i.e., 
less than 5 percent) for both profiles. Likewise, 
the values of vacuum specific impulse, c*, and 
Cp were all within 0.17 percent of the baseline 
calculations for both wall temperature profiles. 
The displacement thickness predicted by the inte- 
grated procedure using the +500 °R wall tempera- 
ture profile is 14.8 percent greater and the 
-500 R wall temperature profile is 15.0 percent 
less than the original predictions (Fig. 12). 

The effect of the +500 R temperature profile on 
heat loss rate resulted in a 7.8 percent exit 
plane value decrease from the baseline value and 
the -500R wall temperature profile resulted in a 
6.1 percent exit plane value increase (Fiq. 13). 

Based on this comparison, the prediction by 
the integrated procedure of the momentum thick- 
ness, local skin friction coefficient, and thrust 
decrement is seen to be relatively insensitive to 
wall temperature profile, whereas the displacement 
thickness and heat loss predictions are very sen- 
sitive. This indicates that if the analysis being 
done is for nozzle design, the wall temperature 
is a very critical input. A nozzle with too great 
or too small of a potential wall displacement will 
result in different exit pressures, Mach numbers, 
and wall angles than those desired. However, if 
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the analysis is to determine specific impulse or 
skin friction of an already fabricated nozzle, an 
approximate temperature profile can be used to 
yield satisfactory results. 

Concluding Remarks 

An analytical study was conducted which com- 
pared the outputs from two different procedures, 
both of which followed the standard JANNAF meth- 
odology of predicting liquid rocket engine 
performance. The predicted values by both proce- 
dures of characteristic velocity and thrust 
coefficient are within 0.6 and 3.0 percent of 
each other, respectively. Predicted vacuum speci- 
fic impulse values are within 1.3 percent of each 
other. Comparing these parameters alone, the two 
procedures displayed good agreement between 
themselves. 

However, the procedures’ predictions of dis- 
placement thickness and skin friction coefficient 
were within a maximum difference of approximately 
10 percent, momentum thicknesses were within 
20 percent, heat loss rates were within 
30 percent, and thrust loss due to boundary layer 
growth was within 40 percent. This was due to 
the differences between the two procedures' input 
wall pressure profiles and physical models (e.g, 
Cebeci-Smith eddy-viscosity model). Upon inspec- 
tion of these parameters, agreement between the 
two procedures can be seen to be poor. Experi- 
mental test results are required to pinpoint 
which (if either) of the two procedures better 
models the real situation. 

Therefore, if approximate values of the gross 
performance parameters are the object of the anal- 
ysis, either procedure will produce a value within 
current measurement accuracy limits of the other. 
It can also be seen, then, that if accurate axial 
profiles of the boundary layer quantities are 
desired, the procedure used (or even the choice 
of boundary layer modules) will greatly affect 
the result. The disjoint procedure predicts 
greater values of displacement thickness, and 
momentum thickness, and lower values of skin fric- 
tion coefficient, thrust decrement, and heat loss 
rate profiles, as compared to the integrated 
procedure. 

Effect of wall temperature profile on repor- 
ted parameters was investigated by displacing the 
baseline wall temperature profile by ±500 *R. 

Less than ±0.2 percent change was experienced by 
values of characteristic velocity, thrust coeffi- 
cient, vacuum specific impulse, and skin friction 
coefficient. Parameters which showed less than 
±5 percent change included momentum thickness and 
thrust decrement. Greatly affected by the change 
in wall temperature profile were boundary layer 
displacement thickness (changed by • 15 percent) 
and heat loss rate (which displayed a ±7 percent 
change) . 

Other areas of concern deal with the opera- 
tion of the computer programs themselves. Com- 


pletion of an analysis of one rocket nozzle for 
the disjoint procedure required considerably more 
clock time than did one run for the integrated 
procedure. This was due to the manual transmis- 
sion and reconditioning of data from the output 
of one computer program to the input of the next 
computer program. Moving the data between pro- 
grams in the disjoint procedure was done as care- 
fully as possible and it was assumed that no 
error was introduced because of incorrect data 
entry. However, this resulted in increasing the 
complexity of the analysis and prolonged the 
amount of time required to do an entire perform- 
ance prediction analysis. Automated interfacinq 
between programs was a major step forward with 
the introduction of TDK. 85 (i.e., the integrated 
procedure). The integrated process eliminated 
this possible source of error and also reduced 
the total time required for a thorouqh analysis. 
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APPENDIX A 


INVISCID CONTOURS 


E = 

60:1 

E = 200:1 

E = 

400:1 

E = 1000:1 
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z. 

R. 

z. 

R, 

Z, 
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in. 

in. 

in. 

in. 

in. 

in. 

in. 

0.50000 

0.00000 

0.50000 

0.00000 

0.50000 

0.00000 

0.50000 

0.00000 

.52980 

.10597 

.55153 

.13874 

.55388 

.13949 

.55639 

.14022 

.56325 

.16766 

.59457 

.20374 

.60268 

.20635 

.61461 

.20968 

.59797 

.22951 

.63923 

.26926 

.65319 

.27376 

.67446 

.27976 

.63416 

.29209 

.68588 

.33600 

.70587 

.34248 

.73637 

.35120 

.67215 

.35607 

.73490 

.40467 

.76129 

.41340 

.80116 

.42507 

.71218 

.42222 

.78640 

.47662 

.81966 

.48771 

.87004 

.50288 

.75426 

.49182 

.84161 

.55344 

.88227 

.56736 

.94477 

.58677 

.79958 

.56659 

.90207 

.63747 

.95098 

.65477 

1.0277 

.67931 

.84938 

.64881 

.96987 

.73163 

1.0283 

.75311 

1.1196 

.78382 

.90556 

.74142 

1.0459 

.83929 

1.1155 

.86613 

1.2243 

.90482 

.96802 

.84720 

1.1331 

.96480 

1.2160 

.99859 

1.3461 

1.0477 

1.0390 

.97000 

1.2343 

1.1129 

1.3331 

1.1558 

1.4883 

1.2185 

1.1210 

1.1142 

1.3497 

1.2880 

1.4679 

1.3429 

1.6536 

1.4235 

1.2121 

1.2820 

1.4821 

1.4949 

1.6232 

1.5652 

1.8453 

1.6690 

1.3150 

1.4776 

1.6303 

1.7361 

1.7988 

1.8263 

2.0638 

1.9598 

1.4255 

1.6998 

1.7915 

2.0114 

1.9917 

2.1265 

2.3063 

2.2972 

1.5447 

1.9494 

1.9648 ! 

2.3201 

2.1993 

2.4644 

2.5693 

2.6798 

1.6673 

2.2218 

2.1468 

2.6604 

2.4194 

2.8396 

2.8505 

3.1077 

1.7947 

2.5208 

2.3386 

3.0379 

2.6539 

3.2588 

3.1532 

3.5903 

2.0089 

3.0658 

2.6997 

3.8075 

3.1041 

4.1248 

3.7441 

4.6011 

2.1994 

3.6005 

3.0833 

4.7110 

3.5887 

5.1512 

4.3913 

5.8145 

2.4603 

4.4208 

3.4157 

5.5728 

4.0157 

6.1401 

4.9708 

6.9967 

2.6743 

5.1883 

3.8634 

6.8626 

4.6000 

7.6336 

5.7804 

8.8049 

2.9330 

6.2492 

4.3107 

8.3208 

5.1980 

9.3425 

6.6255 

10.898 

3.1450 

7.2652 

4.6689 

9.6290 

5.6854 

10.888 

7.3290 

12.812 

3.3071 

8.1523 

5.0843 

11.333 

6.2638 

12.919 

8.1803 

15.351 

3.5316 

9.6024 

5.4424 

12.994 

6.7733 

14.915 

8.9475 

17.870 

3.7040 

10.965 

5.9134 

15.518 

7.4628 

17.977 

10.013 

21.772 

3.8729 

12.641 

6.4260 

6.8083 

7.0711 

18.868 

21.976 

24.576 

8.2451 

8.8592 

9.3302 

22.086 

25.944 

29.418 

11.267 

12.295 

13.116 

14.326 

15.811 

27.075 

32.116 

36.702 

44.635 

57.196 
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TABLE 1. - INPUTS REQUIRED FOR THE DISJOINT PROCEDURE 


Program 

Inputs 

From 

ODE 

Chamber pressure 
Mixture ratio 

Area ratio ' 

Fuel and oxydlzer enthalpies 

Spec 

If led 

OOK 

Mole fractions at throat 
Throat temperature 
Throat velocity 
Pressure at throat 

01 

i 

)E 

BLIMP 

Contour points 
Throat area 
Chamber pressure 

Wall temperature profile 

Spec 

fled 


1 wall pressure profile 

ODK 


! Wall contour points 

Specified 


Throat radius 

Specified 


Chamber pressure 

Options used In BLIMP program: 

Nodal refit on. 

Cebecl-Smlth turbulence model. 
Buddenberg-WIlke mixture formula 
for viscosity. 

Mason-Saxena with Eucken correction 
for thermal conductivity. 

Specified 


Table 2. - INPUTS REQUIRED FOR THE INTEGRATED PROCEDURE 


Program 

Inputs 

From 

TDK. 83 

Contraction ratio 
Radius of throat 

Upstream throat radius of curvature 

Downstream throat radius of curvature 

Wall angle of convergent section 

Radius of curvature Into convergent section 

Nozzle exit angle 

End-of-clrcle angle 

Wall contour points 

Area ratio 

Chamber pressure 

Mixture ratio 

Specified 
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TABLE 3. - LIST OF CHEMICAL 
REACTIONS AND THREE-BODY 
RECOMBINATION-DISSOCIATION 


REACTIONS 


(a) Chemical reactions 



dissociation 

reactions 


H + 0 2 = H 2 0 + M 

H + H = H 2 + M 

M + H 2 0 2 = 20H 

H + OH = H 2 0 + M 


TABLE 4. - NOZZLE EXIT PLANE 
RADII AND LENGTHS 


Area ratio 

60:1 

200:1 

400:1 

1000:1 

Radius (In,) 
Length (In.) 

3.87 

12.64 

7.07 

24.58 

10.00 

35.51 

15.81 

57.20 


TABLE 5. - LIST OF ROCKET ENGINE GEOMETRY PARAMETERS 


ECRAT 


Contraction ratio of combustion chamber 

= 4.223 

RW1U 

- 

Upstream throat radius of curvature 

= 2.0 

RW1D 

_ 

Downstream throat radius of curvature 

= 0.4 

RI 

- 

Normalized Inlet wall radius 

= 1 E-4 

EPS 

- 

Nozzle expansion ratio (see table for values) 


THETA 


Nozzle attachment angle (see table for values) 


THETAI 

- 

Nozzle Inlet angle (see table for values) 


THE 


Nozzle exit angle (see table for values) 



Area ratio 


Angle 

60:1 

200:1 

400:1 

1000:1 

THETAI 

25.00 

25.00 

25.00 

25.00 

THE 

5.05 

5.26 

5.47 

5.74 

THETA 

27.73 

32.61 

35.38 

39.41 
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TABLE 6. - NUMBER OF POINTS 
GIVEN TO SPECIFY CONTOURS 


Procedure 

60:1 

200:1 

400:1 

1000:1 

BLIMP 

38 

40 

40 

40 

TDK 

43 

8 

8 

11 

Initial line points/number of segments 

TDK 

100/3 

100/5 

200/5 

200/5 


TABLE 7. - EXIT PLANE VALUES OF BOUNDARY LAYER PARAMETERS 


[Entries are disjoint/integrated procedure values. Numbers In parentheses are 
(+500 °R/-500 °R) wall temperature profile values using Integrated procedure.] 


Area ratios 

Parameter 

60:1 

200:1 

400:1 

1000:1 

Displacement 
thickness (ft) 

0.00448/0.00455 

0.01755/0.01587 
( .01 822/. 01348) 

0.03560/0.02551 

0.08701/0.06409 

Momentum 

thickness (ft) 

.00430/. 00385 

.00851/. 00764 
(.00737/. 00774) 

.00851/. 01054 

.02045/. 01760 

Thrust 

decrement ( 1 bf ) 

36.76/42.03 

41 .21/51 .74 
(49.38/53.04) 

43.29/55.74 

45.43/66.96 

Skin friction 
coefficient 

.00316/. 00338 

.00294/. 00304 
(.00303/. 00304) 

.00279/. 00311 

.00270/. 00274 

Heat loss rate 
(BTU/sec) 

592.6/696.0 

631.4/789.0 

(727.7/837.0) 

655.1/847.5 

678.7/927.4 
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TABLE 8. - PERFORMANCE COMPARISON 


[Constants: g = 32.1739 ft/sec/sec; 
pc - 1000 psla; At = 0.00545415 ft?.] 


(a) The disjoint procedure 


Area 

ratio 

C*(0DE) 

ft/sec 

C*( BL) 
ft/sec 

C F (OOE) 

Cf(BL) 

60:1 

200:1 

400:1 

1000:1 

7575 

7576.41 

1 .877 
1 .977 
2.020 
2.063 

1 .7852 
1.8577 
1 .8934 
1 .9221 


(b) The Integrated procedure 


Area 

ratio 

C*(0DE) 

ft/sec 

C*(BL) 

ft/sec 

C F (0DE) 

C f (BL) 

C F (BL) 

+ 500 °R 

-500 °R 

60:1 

7590 

7622.57 

1 .878 

1 .8723 



200:1 





1 .978 

1 .9313 

1 .9299 

1 .9346 

400:1 





2.021 

1 .9561 



1000:1 





2.065 

1 .9694 




(c) The disjoint procedure 


Area 

ratio 

I s p( ODE ) , 
sec 

dT, 

Ibf 

dI sp 

sec 

I S p(BL), 

sec 

60:1 

459.8 

36.76 

11.02 

444.4 

200:1 

477.7 

41.21 

12.35 

459.9 

400:1 

485.3 

43.29 

12.98 

466.4 

1000:1 

492.9 

45.43 

13.62 

472.7 


(d) The Integrated procedure 


Area 

ratio 

Isp(ODE) 

(sec) 

dT 

(Ibf) 

dlsp 

(sec) 

I S p(BL) 

(sec) 

I S p(BL) 

+500 °R 

-500 °R 

60:1 

461.0 

42.026 

12.68 

443.95 



200:1 

479.0 

51.735 

15.60 

457.55 

457.06 

458.16 

400:1 

486.5 

55.745 

16.82 

463.44 



1000:1 

494.2 

66.959 

20.20 

466.60 
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INPUTS 


PROGRAM 


OUTPUTS 


CHAMBER PRESSURE 
MIXTURE RATIO 
AREA RATIO 
FUEL AND OXIDIZER 
ENTHALPIES 


THROAT AREA 
NOZZLE WALL CONTOUR 


WALL TEMPERATURE PROFILE 


ODE 

H 



f 

ODK 



/THROAT MOLE FRACTIONS 
T THROAT TEMPERATURE 
THROAT VELOCITY 
THROAT PRESSURE 

I 


WALL PRESSURE PROF I LI 
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SPECIFIC IMPULSE 
THRUST COEFFICIENT 
CHARACTERISTIC VELOCITY/ 
DISPLACEMENT THICKNESS 
MOMENTUM THICKNESS 
HEAT LOSS RATE 
LOCAL SKIN FRICTION 
COEFFICIENT 


FIGURE 1. - DISJOINT PROCEDURE FLOW DIAGRAM. 


0 = THETA 
0j = THETA I 



FIGURE 2. - VARIABLES WHICH DESCRIBE NOZZLE GEOMETRY. 



FIGURE 3. - BASELINE WALL TEMPERATURE PROFILE FOR 
1000:1 AREA RATIO NOZZLE. 



FIGURE A. - COMPARISON OF BASELINE, RL10, AND SSME 
WALL TEMPERATURE PROFILES FOR THE 200:1 NOZZLE. 
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LOCAL SKIN FRICTION 

COEFFICIENT/ Cf DISPLACEMENT THICKNESS, FT 





0 1 2 3 M 5 

AXIAL DISTANCE, FT 


FIGURE 7. - COMPARISON OF LOCAL SKIN FRICTION COEFFICIENT 
NOZZLE AREA RATIO = 1000:1. 



0 12 3 4 5 


AXIAL DISTANCE, FT 

FIGURE 8. - COMPARISON OF THRUST DECREMENT DUE TO BOUND- 
ARY LAYER GROWTH, NOZZLE AREA RATIO = 1000:1. 
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HEAT LOSS# Btu/sec 


PROCEDURE 


1000 

800 

600 

900 

200 



500 



INTEGRATED# EQUILIBRIUM I sp 
DISJOINT# EQUILIBRIUM# I 
DISJOINT# BL I sp PERFORMANCE 
INTEGRATED# BL I sp PERFORMANCE 



FIGURE 9. - COMPARISON OF HEAT LOSS TO THE WALL, NOZZLE 
AREA RATIO = 1000:1. 


FIGURE 10. - COMPARISON OF INTEGRATED AND DISJOINT PRO- 
CEDURES I sp PREDICTIONS. 



FIGURE 11. - COMPARISON OF BASELINE WALL TEMPERATURE 
PROFILE# BASELINE + 500 °R# AND BASELINE - 500 °R. 
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0 .5 1.0 1.5 2.0 

AXIAL DISTANCE/ FT 


FIGURE 12. - COMPARISON OF BOUNDARY LAYER DISPLACE- 
MENT THICKNESS BETWEEN BASELINE/ + 500 °R, AND 
- 500 °R WALL TEMPERATURE PROFILES. 



0 .M .8 1.2 1.0 2.0 

AXIAL POSITION, FT 

FIGURE 13. - COMPARISON OF HEAT LOSS RATE TO THE WALL 
BETWEEN BASELINE, + 500 °R, AND - 500 °R WALL TEMPERA- 
TURE PROFILES. 
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